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Abstract 

We make an analysis of the hadronic invariant mass spectrum in inclusive 
charmless semileptonic B meson decays in a QCD-based approach. The decay 
width is studied as a function of the invariant mass cut. We examine their 
sensitivities to the parameters of the theory. The theoretical uncertainties 
in the determination of |V^| from the hadronic invariant mass spectrum are 
investigated. A strategy for improving the theoretical accuracy in the value 
of \Vub\ is described. 
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I. INTRODUCTION 



The Cabibbo-Kobayashi-Maskawa (CKM) matrix element V U b is a fundemental parame- 
ter in the standard model of particle physics, which is not yet known with sufficient precision. 
A precise determination of it is very important for testing the standard model and for study- 
ing the origin of CP violation and quark masses. 

Until recently, the sole evidence for charmless semileptonic B meson decays and for a 
non-zero value of V u b, which occurs as a multiplicative factor in the decay amplitude, was 
reported from the inclusive charged lepton (electron or muon) energy spectra [|l],|2| . Since the 
rate of the charmless semileptonic B decay is very small, the main experimental challenge 
in these analyses is the suppression of backgrounds. The charged lepton from the inclusive 
charmless semileptonic decay B — > X u iv (£ = e or /i) can be more energetic than that from 
the dominant semileptonic b —>■ c decay; for a B meson at rest, the spectrum in the endpoint 
region (M/2)(l - M 2 D /M 2 ) < E e < (M/2)(l - M 2 /M 2 ) (M being the mass of the B meson, 
Md the mass of the D meson, and M n the mass of the pion) results only from the b — > u 
decay. This simple kinematic cut has been used to obtain the first sign of semileptonic b — > u 
transitions^ and determine \V u b\ from the inclusive charged lepton energy spectra H. The 
theoretical description of the decays requires solving QCD in the nonperturbative regime. 
The accuracy of this inclusive determination of \V u b\ is limited by model dependence, but 
recent work with a QCD-based approach shows promise of improving the situation ||. 

Recently, the CLEO Collaboration has studied the exclusive charmless semileptonic B 
decays to niu, p£u, and loiv |J. This is achieved through the reconstruction of semileptonic 
decays involving a b —>■ u transition. The final states are reconstructed using an observed 
hadron, an e or /i, and a neutrino. Using the hermeticity of the CLEO detector, they re- 
construct the neutrino by inferring its four-momentum from the missing energy and missing 
momentum in each event. With this neutrino reconstruction technique, they have measured 
the branching fractions for B° — > tt~£ + i> and B° — > p~£ + v and have extracted \V u b\ from 
this study. Current determination of \V u b\ by this exclusive method has large model depen- 
dence One should notice, however, recent significant theoretical progress in describing 
exclusive semileptonic b — > u decays 

Various other methods to determine | V u b\ have been proposed. These efforts are obviously 
necessary in order to find the most precise way to extract this important parameter. At 
least, different methods and analyses including two existing types of measurements lend 
confidence and indicate limitations. Among them it has been proposed || long time ago 
that the hadronic invariant mass spectrum in the inclusive semileptonic B decays may be 
useful for the determination of \V u b\. For a B meson at rest, events with a kinematic cut 
on the hadronic invariant mass in the final states, Mx < Mp, come only from the b — > u 
transition (the full b — > u kinematic range is M n < Mx < M) and are thus separated 
from those due to the b — > c transition. Using this cut we have access to a much larger 
region of phase space than with the cut on E e . About 90% of the b — ► u events survive 



1 A new inclusive analysis was presented by ALEPH which is the first evidence for semileptonic 
b — > u transitions in 6-hadrons produced at LEP Q . 
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this cut in contrast to the charged lepton energy cut, with which only about 10% of the 
b — > u events survive. Hence the theoretical calculation of the decay rate is expected to 
be more reliable with the hadronic invariant mass cut. Recent experimental work on the 
neutrino reconstruction technique, which has been used by CLEO in the exclusive charmless 
semileptonic B decay measurement 0, shows promise of measuring the hadronic invariant 
mass spectrum. A number of recent papers [10 have renewed discussions on such possibility 
of determining \ V u b\. On the other hand, studies of the hadronic invariant mass spectra will 
also provide a testing ground for our understanding of the decay dynamics. 

In this paper we will study the hadronic invariant mass spectrum and the decay width 
with a hadronic invariant mass cut in inclusive charmless semileptonic B meson decays in 
a QCD-based approach jll[l2||. This approach describes inclusive semileptonic B meson 
decays by making use of the light-cone expansion and the heavy quark effective theory 
(HQET) and has been recently extended to inclusive semileptonic decays of an unpolarized 
6-flavored hadron |T3]| . The light-cone dominance attributes the nonperturbative QCD 



effects on the underlying weak decays to a single distribution function. The sum rules derived 
from operator product expansion and heavy quark effective theory constrain the distribution 
function considerably. The interplay between nonperturbative and perturbative QCD effects 
has been accounted for and a coherent treatment of both has been formulated. This approach 
has been tested on inclusive semileptonic decays of B mesons. The resulting lepton energy 
spectrum is in good agreement with the experimental data ||. The calculation of the 
semileptonic decay width of the B meson has shown |12| that the kinematically enhanced 
nonperturbative QCD contributions play numerically an important role and \V c b\ has then 
been determined, which is consistent with the exclusive measurements ||14|| . It has been 



shown [0 that an improved precision in \V u b\ can be gained from the lepton energy spectra 
with a controlled theoretical uncertainty. In this work we turn our attention to investigate 
the theoretical accuracy of the alternative determination of \V u b\, i.e., from the hadronic 
invariant mass spectrum. 

Section II presents the theoretical methods to calculate the charmless hadronic invariant 
mass spectrum. We include both the bound-state effects and the QCD radiative corrections. 
In Sec. Ill we examine the sensitivities of the spectrum to the parameters of the theory. The 
effect of the QCD radiative corrections is also discussed. The theoretical uncertainties in the 
calculation of the semileptonic decay width to a given maximum hadronic invariant mass for 
B — > X u £u are studied in Sec. IV. We estimate the theoretical errors in the determination 
of | V u f) | from the hadronic invariant mass spectrum. Finally, concluding remarks and a 
discussion are given in Sec. V. 



II. THEORETICAL METHOD 

The approach has been described in [11-13,5], which we refer to for a more complete 
exposition of the theory. We begin by briefly reviewing the methods relevant to this analysis. 
Because of the light-cone dominance, the nonperturbative QCD effects on the inclusive 
semileptonic B meson decays can be disentangled and encoded in a distribution function: 

B) |, 2=0 , (2.1) 
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where P denotes the four-momentum of the B meson. Several important properties of the 
distribution function were derived from QCD. Due to current conservation, it is exactly 
normalized to unity with a support < £ < 1: 



Jo 



(2.2) 



It obeys positivity. It contains the free quark decay as a limiting case with /(£) = <5(£ — 
rrib/M), where m& labels the b quark mass. In addition, two sum rules for the distribution 
function can be obtained by using the operator product expansion and the HQET method 



15]. These two sum rules determine the mean value ji and the variance a of the distribution 



function, which characterize the position of the maximum and the width of it, respectively: 
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with Eb = G b + K b . The first matrix element G b measures the chromomagnetic energy due 
to the spin coupling between the b quark and the light constituents in the B meson and is 
determined by the mass splitting between B* and B mesons. For the observed difference 
Mb* -M b = 0.046 GeV ||], one gets 



m b G b 



M, 



B 



-0.034 GeV. 



(2.7) 



The second matrix element Kb measures the kinetic energy of the b quark in the B meson. 
A calculation using QCD sum rules yields fL? 



2m\K h = 0.5 ± 0.2 GeV 2 . 



(2.8) 



Taking m& = 4.9 ± 0.2 GeV, the mean value and the variance of the distribution function 
are estimated to be 



fj. = 0.93 ±0.04 



a 2 = 0.006 ±0.002. 



(2.9) 
(2.10) 
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A non-zero value for a 2 indicates the deviation from /(£) = 5(£ — m&/M), i.e., the free quark 
decay, and is a measure of the amount of the nonperturbative QCD contribution. 

To obtain the hadronic invariant mass distribution dT/dMx, it is technically convenient 
to start from the double differential decay rate d 2 T /dq 2 dq° , since the momentum transfer to 
the lepton pair, q, is a quantity which has no distinction between the hadron level and the 
quark level. Following ||, one can account for the nonperturbative and perturbative QCD 
effects on d 2 T/dq 2 dq° in a coherent way through the convolution of the quark- level decay rate 
and the distribution function. This double differential decay rate is then straightforwardly 
converted to dT/dMx- In the B meson rest frame, the hadronic invariant mass spectrum in 
the inclusive charmless semileptonic B decay is finally given by 
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(2.11) 



where p& denotes the four-momentum of the b quark. The detailed expressions of the quark- 
level double differential decay rate d 2 Tf,/dq 2 dq° including the perturbative QCD corrections 
can be found in [18|. The masses of the leptons and the u quark have been neglected. To 
illustrate the effect of the radiative QCD corrections, it is instructive to retain only the terms 
which do not arise from such gluon radiation. Subtracting the radiative QCD corrections, 



Eq. (|2.11|) then reduces to 
dT 
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By integrating Eq. Q2.11|) over the hadronic invariant mass Mx, we obtain the decay 
width for the B — > X u iv decay to a given maximum hadronic invariant mass M™' defined 
by 



r(M; 
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The known properties of the distribution function stated above improve the theoretical 
accuracy remarkably. However, since the distribution function has not yet been completely 
determined in QCD, for pratical calculations we shall adopt the following parametrization 



1 12]] of the distribution function which incorporates all known properties: 
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(2.15) 



where a, b, a, and (5 are four parameters, which are constrained by the sum rules ( |2.3|) and 
(|2.4|) , and iV is the normalization constant. For a — (3 — 1, Eq. (|2.15|) becomes the ansatz 
of [[□]]. In the following, a — /3 — 1 is preset unless explicitly statedQ. 



2 In general, the parameters a and j3 need not be integer. 
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FIG. 1. The charmless hadronic invariant mass spectrum with (solid curve: a s = 0.25, dashed 
curve: a s = 0.30) and without (dotted curve) radiative QCD corrections. The mean value and the 
variance of the distribution function are kept fixed to be /j, = 0.93 and a 2 = 0.006. 

III. CHARMLESS HADRONIC INVARIANT MASS SPECTRUM 

We proceed to examine the sensitivity of the charmless hadronic invariant mass spectrum 
to the parameters of the theory. We calculate the hadronic invariant mass spectrum in the 
inclusive charmless semileptonic B meson decay for a B meson at rest by use of Eqs. ( |2.11|) 



and Q2.15| ). The input parameters are the strong coupling constant a s and the parameters 
which arise in the distribution function. 

The effects of the radiative QCD corrections to the spectrum are demonstrated in Fig. [l|. 
We calculate the spectra using two different values of the strong coupling constant: a s = 0.25 
(solid curve) and a s = 0.30 (dashed curve). The spectrum appears to be insensitive to the 
value of a s , varied within a reasonable range. For comparison, the spectrum without the 
radiative QCD corrections (dotted curve) is also shown, calculated by use of Eqs. ( |2.12| ) 



and ( 2.15 ). It is evident that the charmless hadronic invariant mass spectrum receives large 
radiative QCD corrections. This is easily understood because gluon bremsstrahlung strongly 
affects the hadronic final states. 

We explore next the sensitivity of the spectrum to the form of the distribution function. 
The mean value and the variance of it vary in the ranges specified in Eqs. ( [2.9|) and ([2.1 0|) 



inferred from the sum rules (2J3) and (|2.4| ). We show in Fig. ^| the variation of the hadronic 
invariant mass spectrum due to the changes in the mean value fi and the variance a 2 of 
the distribution function. The position of the maximum for the spectrum is a sensitive 
function of the mean value /i. The sensitivity of the spectrum to the variance a 2 is seen to 
be relatively weak. 

Furthermore, to get an idea of the sensitivity of the spectrum to the form of the distribu- 
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FIG. 2. Dependence of the charmless hadronic invariant mass spectrum on the mean value and 
the variance of the distribution function. The strong coupling constant is taken to be a s = 0.25. 



tion function when keeping the mean value and the variance of it fixed, we change the shape 
of the distribution function by varying the values of the two additional parameters a and (3 
in the parametrization (|2.15|) (for given values of a and (3, the parameters a and b are fixed 
by the mean value and the variance). In Fig. |3| we plot the distribution functions for four 
different choices of the parameters a and (3 in ( |2.15| ) that give the same mean value /i = 0.93 
and variance a 2 = 0.006. The spectra calculated using the two very different distribution 
functions with a = 1 and (3 = 1 (solid curve), a = 2 and (3 = 2 (dashed curve), respectively, 
at fixed mean value and variance are shown in Fig. |j. It is seen that the hadronic invariant 
mass spectrum relies strongly on the form of the distribution function, even when the mean 
value and the variance of it are kept fixed. 



IV. PARTIAL (TOTAL) DECAY WIDTHS AND DETERMINATION OF \V ub \ 

Using the charmless hadronic invariant mass spectrum, the semileptonic branching frac- 
tion A£> M (M^ ut ) for the charmless semileptonic B decay with the hadronic final states up to a 
given maximum invariant mass M™* can be measured. Together with the measured lifetime 
of the B meson, tb, and the theoretically calculated decay width defined in Eq. (|2.14j ) up 
to the CKM factor IKjJ 2 , this determines \V u b\- 



|2 A0 u (M' 



cut^ 
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cut 



where 7„(M^ ut ) is defined by T(M^ ut ) = |T4b| 2 7 u (^ ut ) 
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FIG. 3. The distribution functions from the parametrization ( |2.15| ) for four different pairs of 
the parameters a and (3. All the distribution functions give fi = 0.93 and a 2 = 0.006. 
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FIG. 4. Comparison of the charmless hadronic invariant mass spectra in two distribution func- 
tions. The mean value and the variance of the distribution functions are fixed to be fi = 0.93 and 
a 2 = 0.006. We take a s = 0.25. 
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We investigate the theoretical uncertainties in this determination of \V u b\. The semilep- 
tonic decay width T(M™ t ) is displayed in Fig. [|as a function of the hadronic invariant mass 
cutoff, M£ ut . For M^ ut = M, r(M£ ut ) is the total semileptonic decay width for B -> XJu. 
The solid and dotted curves correspond to two different values of a s : a s = 0.25 (solid), 
a s = 0.30 (dotted), and the other parameters are identical. We find that a reasonable 
variation of the value of a s leads to a small change in the value of the decay width r(M^ ut ). 

In order to explore the impact of the form of the distribution function, as in the last 
section, we first study the sensitivities of r(M^ ut ) to the mean value and the variance of 
the distribution function by using the parametrization ( |2 . 1 5 ) setting a = (3 = 1. We vary 
the mean value and the variance in the ranges of Eqs. ( |2.9| ) and ( |2.10| ). The resulting decay 
widths are shown in Fig. |] corresponding to /i = 0.93 and cr 2 = 0.006 (solid curve), \i = 0.89 
and a 2 = 0.006 (long-dashed curve), /i = 0.93 and cr 2 = 0.008 (dot-dashed curve), with 
a s = 0.25 held fixed. Then we change the form of the distribution function by changing 
the values of the two additional parameters a and f3 in Eq. ( [2.15| ) in order to examine the 
further sensitivity of T(M™ t ) to the distribution function when keeping the mean value 
and the variance of it fixed. The results are also shown in Fig. |5| for a — (3 — 1 (solid 
curve), a = P = 2 (short-dashed curve), with \i = 0.93, cr 2 = 0.006, and a s = 0.25 held 
fixed for both curves. We find that for the cut M^ ut less than about 1.2 GeV, the partial 
decay width r(M^ ut ) is rather sensitive to the form of the distribution function, whereas for 
M^ ut greater than about 1.2 GeV, the decay width r(M™ t ) (including the total semileptonic 
decay width for B — > X u iv decays) is sensitive essentially only to the mean value of the 
distribution function, which is calculable from HQET with the results given in Eqs. ( |2.3| ) 
and ( |2.9| ). These imply that the theoretical uncertainties in the determination of \V u b\ are 
under control for the value of the hadronic invariant mass cutoff greater than about 1.2 GeV, 
but very large for the value of the cutoff less than 1.2 GeV or so. 

Adding the errors on the semileptonic decay width due to the value of a s and the form 
of the distribution function linearly to be conservative, we estimate the theoretical errors in 
the determination of \V u b\. Without any cut, the theoretical error in determining \V u b\ from 
the total charmless semileptonic decay width is about 13%. Applying a cut on the hadronic 
invariant mass, the theoretical errors in determining \ V u b\ from the hadronic invariant mass 
spectrum are approximately 16% for M^ ut = 1.8 GeV, 26% for M^ ut = 1.5 GeV, and 44% 
for Mf 11 = 1.2 GeV, respectively. 



V. CONCLUSION AND DISCUSSION 



We have presented an analysis of the hadronic invariant mass spectrum and studied the 
decay width as a function of the invariant mass cut for inclusive charmless semileptonic 
B decays. An approach based on the light-cone expansion and the heavy quark effective 
theory has been exploited. Nonperturbative QCD effects in the processes are contained in a 
distribution function and several important properties of it are known from QCD. However, 
one needs a parametrization of the distribution function incorporating the known properties 
in order to make quantitative predictions, since it cannot yet been completely determined 
from QCD. For our analysis, we employ the distribution function parametrization given in 
Eq. (|2.15| ). We have examined the sensitivities of the hadronic invariant mass spectrum and 
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FIG. 5. The semileptonic decay width defined in Eq. ( 2.14 ) for B — > X u lv as a function of 
the hadronic invariant mass cutoff. The dependence on the mean value and the variance of the 
distribution function is shown using a s = 0.25 and the parametrization ( 2.15| ) with a = j3 = 1 and 



three pairs of the mean value and the variance: \x = 0.93 and a 2 = 0.006 (solid curve), \i = 0.89 
and a 2 = 0.006 (long-dashed curve), \i = 0.93 and a 2 = 0.008 (dot-dashed curve). The dependence 
on the form of the distribution function at fixed mean value and variance is shown by the solid and 
short-dashed curves corresponding to the same values of a s and of the mean value and the variance, 
but a = (3 = 2 for the short-dashed curve. The dependence on the strong coupling constant is 
illustrated by the solid and dotted curves obtained from the same parameters except a s = 0.30 for 
the dotted curve. 



the decay width with a hadronic invariant mass cut to the parameters of the theory, namely 
the strong coupling constant and the form of the distribution function. 

It is shown that a reasonable variation of the value of a s hardly affects the charmless 
hadronic invariant mass spectrum. The dependence of the decay width with a hadronic 
invariant mass cut on the value of a s is small. The main theoretical uncertainties arise from 
the form of the distribution function. 

We found that the charmless hadronic invariant mass spectrum exhibits a significant 
sensitivity to the form of the distribution function. Thus, if experimentally feasible, the 
measured hadronic invariant mass spectrum will impose strong constraints on the shape of 
the distribution function. 

We have also studied the sensitivity of the semileptonic decay width to a given maximum 
hadronic invariant mass to the form of the distribution function. For the hadronic invariant 
mass cutoff greater than about 1.2 GeV, the value of the decay width is sensitive essentially 
only to the mean value of the distribution function whose value is known from HQET, and 
thereby the result is nearly model independent and the uncertainties in the calculation of the 
decay width are under control. However, the theoretical uncertainties in the partial decay 
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width with a hadronic invariant mass cutoff less than 1.2 GeV are very large, obscuring its 
usefulness for determining \V u b\, since it is rather sensitive to the form of the distribution 
function. 

We have observed the feature that for M™ 1 > the decay width tends to show only 
little dependence on the invariant mass cut (see Fig. 5) and most of the b — > u events lie 
below M x = M D (see also Figs. 1, 2, 4). 

We have estimated the theoretical errors in determining \V u b\ in our approach. At the 
present time, the theoretical error on \V u b\ from the total semileptonic decay width for the 
B — > X u lv decays is about 13%. Applying a hadronic invariant mass cutoff, the theoretical 
errors on \Vub\ are approximately 16% for M™* = 1.8 GeV, 26% for M^ ut = 1.5 GeV, 
and 44% for M™* = 1.2 GeV, respectively. The higher the hadronic invariant mass cutoff 
can be experimentally made to be, the smaller the theoretical error on \V u b\; the smallest 
theoretical error on \ V u b\ would be achieved if the total charmless semileptoinc decay width 
can be measured^. 

Besides the hadronic invariant mass cutoff, the mean value of the distribution function is 
another key quantity for precise determination of \V u b\, which requires as precise knowledge 
of the mean value as possible. Based on this analysis and the previous one ||, a detailed fit to 
the measured charged-lepton energy spectrum and/or the hadronic invariant mass spectrum 
can be used to extract the parameters (especially, the mean value of the distribution function) 
experimentally. This procedure allows to reduce the uncertainties in the determination of 
\V u b\. Future measurements of the distribution function and continuing efforts in calculations 
of hadronic matrix elements that govern B decays should enable to further reduce the 
uncertainties. 
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3 An attempt has been made || by the ALEPH Collaboration to measure the inclusive charmless 
semileptonic branching fraction of B hadrons. 
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